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Abstract 

The flapping equation for a helicopter in for- 
ward flight has coefficients which are periodic in 
time, and this effect complicates the calculation 
of stability. This paper presents a constant 
coefficient approximation which will allow the use 
of all the well known methods for analyzing constant 
coefficient equations. The flapping equation is 
first transformed into the nonrotating coordinate 
frame, where some of the periodic coefficients are 
transformed into constant terms. The constant 
coefficient approximation is then made by using 
time averaged coefficients in the nonrotating frame. 
Stability calculations based on the approximation 
are compared to results from a theory which cor- 
rectly includes all of the periodicity. The com- 
parison indicates that the approximation is reason- 
ably accurate at advance ratios up to 0.5. 

Notation 

a blade lift curve slope 

B tip loss factor 

c blade chord 

I blade flapping inertia 

i /T 

kg flapping spring stiffness 

N number of blades 

R rotor radius 

t time, sec 

V forward velocity 

a angle of attack of hub plane 

j*i flapping of ith blade relative to hub plane 

6 vector of rotor degrees of freedom in non- 

rotating coordinates 
3 0 rotor coning angle 

3 lc rotor tilt forward (longitudinal flapping) 

&i s rotor tilt to left (lateral flapping) 

B 2 rotor differential flapping 

y blade lock number, pacRyl 

X eigenvalue or root, nondimensionalized by ft, 



M rotor advance ratio, V (cos a)/ftR 
v flapping natural frequency of rotating blade 
o real part of eigenvalue 
P air density 

azimuth angle, fit 
ft rotor rotational speed 

imaginary part of eigenvalue 
flapping n atural frequency of stationary 
blade, yC7T 

(„) derivative, d( )/d^ 

O derivative, d 2 ( )/dt^ 2 

For helicopter stability and control studies, 
it is desirable to use as simple a math model as 
possible while retaining reasonable accuracy, both 
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to reduce computation effort and to gain insight 
into system behavior. However, for a helicopter in 
forward flight, the rotor flapping motion is 
described by a differential equation having coeffi- 
cients which are periodic in time (azimuth) . This 
fact complicates the solution of the equation, 
requiring methods which use considerable numerical 
computation and which give* little insight. Thus it 
is desirable to find a differential equation with 
constant coefficients (hence an approximation) 
which adequately represents the forward flight 
flapping dynamics of a helicopter rotor. If such 
an equation is found, all of the well known tech- 
niques for analyzing constant coefficient equations 
may be used. 

The flapping equation may be transformed into 
the nonrotating coordinate frame, as done in 
References 1 and 2, where some of the periodicity 
is transformed into constant terms. This result 
suggests that the use of constant coefficients in 
the nonrotating frame will retain some of the 
periodic system behavior. The constant coefficient 
approximation examined herein is made by using time 
averaged coefficients in the nonrotating frame. A 
comparison is made between the eigenvalues (sta- 
bility) obtained from the approximation and the 
results from a theory which correctly includes all 
of the periodicity. The comparison indicates that 
the approximation is a useful representation of 
helicopter flapping dynamics for both hingeless and 
articulated rotors. This approximation was briefly 
discussed in Reference 1 for one set of rotor 
parameters. The present paper discusses the 
approximation in a more general manner and gives 
more insight into its features, limits, and 
applicability. 

The rotor math model used here is for fixed 
shaft operation and includes only first mode 
(rigid blade) flapping, with spring-restrained 
flapping hinges at the hub center. Flapping 
natural frequency may be matched by selecting the 
spring rate. Thus the only approximations are in 
the use of the aerodynamic terms for rigid blade 
motion. Uniform inflow is used, and for the 
advance ratios considered here (y < 0.5), reverse 
flow effects are not included. 

Equations of Motion 

In this section, the single blade homogeneous 
flapping equation is presented for a rigid, spring- 
restrained, centrally hinged blade. This equation 
is then transformed to a nonrotating coordinate 
frame, using a coordinate transformation which is 
briefly discussed. Insight into the fundamental 
behavior of the rotor is gained by examining the 
hovering (y » 0) eigenvalues of the equation in 
nonrotating coordinates. 
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For the single blade, the homogeneous equation 
of motion is 


where 


+ M | B i + (v2 + V 6 i = ° 


M B = I 84 + y <> fi3 Sin *i 
Mg = y 1 B 3 cos + y 2 -g- B 2 sin 2 i|k 
, “ e 2 k s 

v 2 = 1 + 4 _ = 1 + - 1 - 


( 1 ) 


& 


in 2 


Note that reverse flow has not been included here. 
Although it could be included, it would not signifi- 
cantly affect the results for p < 0.5, since the 
additional terms are fourth order in p. 

By a coordinate transformation of the Fourier 
type, the single blade equation may be written in 
terms of nonrotating coordinates. The transforma- 
tion accounts for the motion of all blades, and the 
number of degrees of freedom is equal to the number 
of blades. For example, with a three-bladed rotor, 
the degrees of freedom are coning (all blades flap- 
ping together), rotor pitching (cosine \p flap- 
ping), and rotor rolling (sine \p flapping). 

Adding a fourth blade adds a differential flapping 
degree of freedom, where blades 1 and 3 flap in one 
direction while blades 2 and 4 flap in the other 
direction. This type of differential motion is a 
degree of freedom with rotors having any even num- 
ber of blades. Adding more blades adds degrees of 
freedom which, in the nonrotating frame, warp the 
plane described by the sine ip and cosine \p 
flapping motion. 

The coordinate for the single blade is B^. 

For a three-bladed rotor, the corresponding 
nonrotating coordinates are 



where B 0 , 3, c , and B ls are rotor coning, pitch- 
ing, and rolling motions. For a four-bladed rotor. 



where B 2 is the differential flapping motion 
discussed above. 

In general, the blade degrees of freedom in 
the transformation are 


, N 

i*l 

2 N 


nc 


i E *i cos n *i 


i= i 
N 


ns 


= n E e i sin n *i 


1*1 

N 


S N = i J V-D 1 . N even only 

2 i=l 

Then the motion of the ith blade is 
K 

B i ■ 6 o + E (8„c cos n *i + B ns sin n *i) + e N ( ' X) ; 


n*l 


j (N - 1), N odd 
i (N - 2), N even 


The equations of motion (that is, eq. (1)) must 
also be converted from a rotating to a nonrotating 
frame by a similar procedure. This process is 
accomplished by operating on the equations with the 
summation operators 

wE ( --- h IE ( --- )C0S n, "i ■ 

i i 

|£(...;)sin n<|/. . Jf^C.-OC-D 1 


This is virtually the same procedure used in 
Reference 1. 

It may be seen that the transformation 
involves multiplication by sin rp, cos sin 2^, 
cos 2ip, etc. This changes some of the periodic 
terms of the equations in the rotating reference 
frame into constants (plus higher harmonics) due to 
products of periodic terms, and vice versa. 

Performing the indicated operations for N * 3 
yields the following equations for a three-bladed 
rotor. 


I B “ 0 u IT 

0 B 1 * ♦ u B J sin 3^ 2-u B 5 cos 3* 

^ B 3 *2 * u yj B 3 cos 3<> IF v s * n 3'('| 

p 2 -fe B 2 sin 3* - u 2 ^ B 2 cos3V> 

i 2 - 1 ♦ p ^ B 3 cos 3* 


|b 5 *u j |b j sin 3* 

-•i** cos 3* - 1 (b 1 * - J u 2 B 2 ) +u!b 3 sin 3^ 


v 2 - 1 - p 1 B 3 cos 3* 


( 2 ) 


Similarly, operating as above with N « 4, the 
equations of a four-bladed rotor are obtained. 
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For three blades, the hovering equations are iden- 
tical, except that the B 2 equation is then absent. 
Hie B 0 and B 2 equations at hover are completely 
uncoupled and are both identical to that of the 
single blade in rotating coordinates. 


B* + j *% + v 2 B ± = 0 


(1), for y = 0 


B + 1 B 4 B„ + v 2 B„ = Ol 
o 8 o o 

b 2 + 1 b 4 b 2 + v z e 2 = o, 


from equation (4) 


(3) 


The three- and four-bladed rotors have similar 
behavior except for the terms which are periodic in 
The periodic terms are 3/rev for the three- 
bladed rotor, but are 2 and 4/rev for the four- 
bladed rotor. 

The main advantage of the transformed equations 
is that it is easier to express the combined rotor 
and airframe motions because the rotor equations 
are now in a nonrotating reference frame and 
include the motions of all blades. Furthermore, 
rotor motions are more intuitively understood, 
since the degrees of freedom are those seen by an 
observer in or beside the helicopter. 

In the nonrotating coordinates of equa- 
tions (2) and (3), the equations are coupled by 
off-diagonal terms. Note however, these are actu- 
ally independent blades (unless some sort of feed- 
back is added) and the coupling is due to the 
coordinate transformation. 

To gain understanding of these degrees of 
freedom, the hovering (v ■ 0) behavior is examined 
next. The hover equations for four blades are 
given below. 


Eigenvalues of these equations are easily calcu- 



called coning and reactionless modes. The reaction- 
less mode is so named because at hover it produces 
no net reaction at the hub. The equations for 
rotor pitching and rolling are. 


lc 


R»‘ 


1 B 4 0 0 



-2 





B + 




-2 

0 


2 

0 




0 

v 2 - 1 
0 


0 

£b“ 

V 2 - 1 
0 


0 

0 

\»2 


B - 0 



and the characteristic equation is then, 

(a 2 ♦ 1 B 4 A + v 2 - l) 2 + (» * i ‘'Y ■ 0 

The eigenvalues for this equation axe shown in 
figure 2. By analogy to a gyro, these modes will 
be called precession (the lower frequency mode) and 
(4) nutation (the higher frequency mode) . The damping, 
-y/16, is the same as for the single blade of equa- 
tion (1) and for the B 0 and B 2 modes discussed 
above. However, the coordinate transformation has 
resulted in the precession mode frequency being G 
lower than the single blade mode frequency in the 
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Periodic Coefficient Solutions 



Figure 2. Hover eigenvalues of precession and nuta- 
tion modes (v * 1.2, y = 8). 


rotating frame. Similarly, the nutation mode fre- 
quency is ft higher than that of the single blade in 
the rotating frame. 


The coning mode (fig. 1) will excite vertical 
motions of the vehicle, while the precession mode 
will excite pitch and roll motions. Thus vehicle 
responses are more intuitively understood by use 
of the nonrotating coordinates. Also, these equa- 
tions may be used to study feedback control systems 
such as rotor tilting or rotor coning feedback, 
which were discussed in Reference 1. Note (from 
eqs. (1) and (2)) that the performance of such sys- 
tems will depend on the number of blades used, 
since the blade motions become coupled by the feed- 
back terms and the coupling will vary with the 
number of blades. 


To compare the various modes with each other 
and with other theories, it is necessary to trans- 
form all eigenvalues into the same reference frame. 
The obvious choice is the rotating coordinates of 
equation (1) , since most other theories are appli- 
cable to this frame. As may be seen by comparing 
figures 1 and 2, the precession and nutation modes 
may be transformed back into the rotating frame by 
adding and subtracting ft respectively. This 
process results in four identical eigenvalues, as 
expected, since the rotor is composed of four iden- 
tical blades, each described in the rotating frame 
by equation (1). As noted above, the frequencies 
of the 6 q and modes do not change 

The equations for p = 0 have been easily 
solved and the nonrotating coordinate system has 
been presented and discussed. In nonhovering 
flight, however, the equations have periodic 
coefficients, which makes the equations more diffi- 
cult to solve, as well as giving the solutions some 
special characteristics. These will be discussed 
in the next section. 


Floquet Theory 


Eigenvalues of equations such as (1) may be 
found with Floquet theory, as for example in 
References 3 and 4. The equation is integrated for 
one period (tJj = 0 ,..., 2tt) for each independent 
initial condition to obtain the state transition 
matrix. The frequency and damping of the system 
modes are then obtained by taking the logarithms of 
the state transition matrix eigenvalues. 

This technique has been applied here to three 
cases, and the results are shown on figure 3 for 



Figure 3. Floquet theory root loci for varying p; 
single blade in rotating coordinates. 

varying p. Note that as p is increased, the 
frequency (w/ft) decreases, while the damping (a/ft) 
remains constant at -y/16 until the frequency 
reaches an integer multiple of 1/2 /rev. As p is 
increased further, the frequency remains constant 
while the damping both decreases (the upper roots) 
and increases (the lower roots) as shown for cases 
a and c. This behavior may be surprising to those 
accustomed to constant coefficient equations, but 
is typical of periodic systems. The nonsymmetry 
about the real axis is analogous to the behavior of 
a constant coefficient equation root locus when the 
locus meets the real axis. At that point, the 
roots separate (no longer complex conjugates) , one 
becoming less stable and the other becoming more 
stable. With periodic coefficient equations, the 
separation can occur at any multiple of 1/2 the 
frequency of the periodicity. Actually, the con- 
stant coefficient equation is a special case of the 
periodic one, where the frequency in the coeffi- 
cients is zero. This behavior may be seen in more 
detail by plotting the eigenvalues versus p, as in 
figure 4, which again shows results from Floquet 
theory. 


48 



L2 r 


i.o h 


at 

a 


.8 r 


a_ 

b. 


a v = u , y-6 

b. v = i.o, y*6 

c. v*i.o t y~\ 2 


C 





Figure 4. Floquet theory variation of frequency and 
damping with p; single blade in rotating coor- 
dinates. 

Figures 3 and 4 have shown the eigenvalues in 
the rotating coordinate system. These may be 
transposed into the nonrotating system to examine 
the behavior of the nonrotating modes. Choosing 
case c (v = 1.0, y * 12) as an example, the root 
locus is plotted on figure 5. The coning mode has 
the same eigenvalues shown in the two previous 
figures. The nutation and precession modes have 
the same damping, but as mentioned before, their 
frequencies are ft higher and lower, respectively, 
than the coning frequency. 

The regions where the frequency remains constant 
while the damping changes, called critical regions, 
may be illustrated by constructing the y - v 
plane as in figure 6 (and discussed in References 3 
and 4). In the 0/rev region, the behavior is like 
that of a constant coefficient equation when the 
root locus meets the real axis; there are two real 
roots, with order p 2 changes in damping. In. the 
j/rev region, the frequency is exactly half of the 
rotational frequency (fig. 3, case b), and the 
damping changes somewhat more rapidly. In the 
1/rev region (fig. 3, case a) the frequency is the 
same as the rotational frequency (ft), and again 
the damping changes are order p 2 . As previously 
noted, damping is constant at -y/16 outside of the 
critical regions. Note that varying v has little 
effect on the boundaries of the 0/rev and -j/rev 
regions, but as v is increased the 1/rev region 
moves upward. 



Figure 5. Floquet theory root loci for varying y; 
three- b laded rotor in nonrotating coordinates 
(v * 1.0, y = 12; case c). 



In this section, the characteristics of the 
periodic coefficient solutions have been discussed 


Figure 6. y - V plane for single blade in rotat- 
ing coordinates based on Floquet theory. 
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Note that the & 2 equation is not coupled to the 
others and is the same as the 6 2 equation for 
hover; hence it yields only the y = 0 roots. 
Therefore the B 2 equation will not be discussed 
further or included in subsequent figures. The 
3 0 equation has only one y-dependent term, 
coupling it to the 6 1S motion. The pitch and 
roll equations are coupled by both damping and 
aerodynamic spring terms. 

Comparison 

As noted earlier, eigenvalues may be compared 
by adding £2 to the precession frequency and sub- 
tracting £2 from the nutation frequency. In 
examining the constant coefficient approximation, 
any differences in eigenvalues will be due to the 
dropped periodicity. That is, all of the roots 
should approximate those obtained by using Floquet 
theory to solve equation (1). Using the comparison 
method mentioned above, the constant coefficient 
approximation is compared to Floquet theory results 
in figures 7, 8, and 9. The frequency scales have 
been expanded to exaggerate the effects of forward 
speed. Each of the three cases is discussed below. 

Case a. v = 1.1, y = 6. (fig. 7) 


Figure 6. Concluded. 

for nonhovering flight. The next sections will 
discuss an approximation which has constant coeffi- 
cients, yet gives some of the behavior of the peri- 
odic coefficient system. 


Constant Coefficient Approximation 


In equation (1) the periodic coefficients are 
all of the speed (y) dependent terms, and a con- 
stant coefficient approximation yields only the 
hover solution. However, in the nonrotating frame 
of equations (2) and (3), these periodic terms 
have been transformed into constants plus higher 
harmonic periodic terms. This result suggests that 
the primary effects of y may be determined by 
using the average values of the coefficients. The 
constant coefficient approximation thus obtained 
for a four-bladed rotor is given in equation (5). 
The corresponding equation for three blades is 
identical, except that the 6 2 motion is absent. 


0 u j2 B 3 0 


0 £ B 4 2 

i * 3 -2 i 


0 

0 

X 




i ( B - . i „*»*) 


o - i {* - i „’»>) 


' - 1 
0 


(5) 


This case corresponds to a hingeless rotor 
similar to the Lockheed XH-51. For this rotor the 
variations with y of frequency and damping are 
small but significant since the 1/rev critical 
region is encountered (see fig. 6). All three 
modes of the approximation agree well with Floquet 
theory at low advance ratios, where the influence 
of the periodic coefficients is small. As the 
advance ratio is further increased, the precession 
mode displays the same type of behavior as the 
Floquet theory results, but the other two modes do 
not. For the precession mode (and the Floquet 
theory), the frequency becomes constant at 1/rev, 
and the damping then has two values as previously 
discussed. It is useful to examine why the con- 
stant coefficient approximation displays periodic 
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Figure 7. Comparison of constant coefficient 
approximation to Floquet theory (v = 1.1, y = 6). 
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Case b. v = 1 . 0 , y = 6. (fig. 8) 

This case corresponds to an articulated rotor 
having relatively heavy blades, such as might be 
used for a high speed helicopter. This case is 
well removed from critical regions, and there are 
no significant changes in the eigenvalues for the 
u range shown. The constant coefficient approxi- 
mation agrees well with results from Floquet theory. 
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Figure 7. Concluded. 

system (critical region) behavior. In this case, 
the precession roots at hover (p = 0) are very near 
the real axis due to the coordinate transformation. 
As p is increased, the precession roots move 
toward the real axis and then split when they reach 
the axis, as usual with constant coefficient sys- 
tems. Thus the damping both increases (the left 
branch) and decreases (the right branch). 


i 1 1 i i 

0 ♦» *2 -3 .4 .5 

H- 

(b) DAMPING 


Figure 8. Comparison of constant coefficient 
approximation to Floquet theory (v = 1.0, y = 6) . 
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Figure 8. Concluded. 
Case c. v = 1.0, y = 12. (fig- 9) 


This case corresponds to a typical articulated 
rotor with blades similar to many aircraft flying 
today. The Floquet theory indicates that the ^-/rev 
region is encountered at y = 0.215. It is seen 
that the nutation mode is a poor approximation. 
Apparently, the constant coefficient approximation 
is not adequate for higher frequency modes if a 
critical region is encountered. The precession and 
coning modes (combined) , however, do display the 
correct type of behavior: the frequency approaches 

^-/rev and the damping both increases (the precession 
mode) and decreases (the coning mode) . In this 
case, the correct behavior is obtained because two 
modes are involved. As may be seen in figure 9(c), 
the two sets of Floquet roots approach each other, 
meet at y/rev, and split (no longer complex con- 
jugates) . This behavior is approximated by the 
coning and precession modes, but in the approxima- 
tion, the roots remain complex conjugates as shown 
in figure 9(c). The frequency of the precession 
mode does not agree well with Floquet results, but 
its damping is increasing; hence it is of less 
interest. The coning mode agrees well with the 
Floquet results, predicting the reduced damping 
very accurately. 


Perturbation Theory 

Equation (1) has also been studied in 
Reference 5, using a perturbation technique known 
as the method of multiple time scales. Analytic 
expressions are derived for the eigenvalues, with 
expressions valid near and within each of the 
critical regions and ones which are valid away from 
the critical regions. These results are very 
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P 
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Figure 9. Comparison of constant coefficient 
approximation to Floquet theory (v *= 1.0, y * 12). 
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Figure 9. Concluded. 


useful; they give additional insight into the behav- 
ior of periodic systems in general and equation (1) 
in particular. A comparison between the Floquet 
results of the present work and the analytic results 
from Reference 5 indicates that the latter are also 
useful quantitatively. An exception is near the 
^/rev region, where the perturbation solution was 
carried only to order p. It should evidently be 
extended to order y 2 as was the rest of the 
perturbation solution. 

Discussion 

Based on the cases described above, it is 
apparent that the constant coefficient approximation 
may be used to calculate rotor eigenvalues at 
advance ratios up to 0.5. A range of rotor param- 
eters (y and v) have been studied which are repre- 
sentative of most conventional helicopters. The 
lower frequency modes agree well with Floquet 
results and display behavior approximating that of 
the Floquet theory critical regions. Therefore, 
there are many cases where the approximation may be 
used instead of more complicated methods. 

The higher frequency modes of the approxima- 
tion, however, do not display the correct behavior. 
Where these modes are important, for example, in 
using high gain feedback, the approximation should 
be used with caution. 

The perturbation theory of Reference 5 is very 
easy to use for rotor stability calculations. 


However, the solutions are for uncoupled blades in 
the rotating coordinate frame. To account for 
inter-blade coupling (as with certain feedback 
schemes) one must either use another technique 
such as that described herein or rederive the solu- 
tions with the coupling included. 

Conclusion 

Transforming the flapping equation of a heli- 
copter rotor in forward flight into the nonrotating 
coordinate frame results in a set of differential 
equations where some of the periodicity due to 
forward flight is transformed into constant terms. 
Using the time-averaged values of these, i.e., 
dropping the remaining periodicity, gives a con- 
stant coefficient approximation which retains some 
of the periodic effects. Comparison between results 
of the approximation and those of Floquet theory 
indicates that the approximation should be accept- 
ably accurate for calculating flapping stability of 
most helicopters for the advance ratios shown 
herein. Use of the nonrotating coordinates has 
given insight into rotor behavior and indicates how 
the vehicle motion would be affected by the rotor 
modes . 

The higher frequency modes of the approximation 
do not agree well with Floquet theory. Where these 
modes are important for example, in using high gain 
feedback control systems, the approximation should 
be used with caution. 
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